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Abstract 

There has been a large amount of work aimed at understanding 
the stability of nonlinear dispersive equations that support solitary 
wave solutions [201 O [H [191 EHl [211 M [351 ES] . Much of this work 
relies on understanding detailed properties of the spectrum of the 
operator obtained by linearizing the flow around the solitary wave. 
These spectral properties, in turn, have important implications for the 
long-time behavior of solutions to the corresponding partial differential 

equation[ini EllI31inilIElEai231l2Sll2Zlllil3Zll331l33]. 

In this paper we consider periodic solutions to equations of Korteweg- 
Devries type. While the stability theory for periodic waves has re- 
ceived much some attention [H ^ IS CE51 H§1 ED 02] the theory is 
much less developed than the analogous theory for solitary wave sta- 
bility, and appears to be mathematically richer. We prove an index 
theorem giving an exact count of the number of unstable eigenval- 
ues of the linearized operator in terms of the number of zeros of the 
derivative of the traveling wave profile together with geometric infor- 
mation about a certain map between the constants of integration of 
the ordinary differential equation and the conserved quantities of the 
partial differential equation. 
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This index can be regarded as a generalization of both the Sturm 
oscillation theorem and the classical stability theory for solitary wave 
solutions for equations of Korteweg-de Vries type. In the case of a 
polynomial nonlinearity this index, together with a related one intro- 
duced earlier by Bronski and Johnson, can be expressed in terms of 
derivatives of period integrals on a Riemann surface. Since these pe- 
riod integrals satisfy a Picard-Fuchs equation these derivatives can be 
expressed in terms of the integrals themselves, leading to an expres- 
sion in terms of various moments of the solution. We conclude with 
some illustrative examples. 



1 Introduction and Preliminaries 

In the area of nonlinear dispersive waves the question of stability is an im- 
portant one, as it determines what states one is likely to observe in practice. 
In this paper we consider periodic traveling wave solutions to equations of 
KdV type: 

u t + u xxx + (f(u)) x = (1) 

where / is assumed to be C 2 . The main goal of this paper is to determine an 
index expressible in terms of the conserved quantities restricted to the man- 
ifold of periodic traveling wave solutions which yields sufficient information 
for the orbital and spectral stability of the underlying wave. Our results are 
geometric in nature as they relate to Jacobians of various maps which relate 
naturally to the structure of the underlying linearized operator, and will be 
most explicit when / is a polynomial. Note that when / is a polynomial 
of degree at most three the above is integrable via the inverse scattering 
transform, but this is not the case for polynomials of higher degree. 

Assuming a traveling wave of the form u = u{x — ct) one is immediately 
led to the following nonlinear oscillator equation 

^- C -^ + au + F(u) = E (2) 

where F is the antiderivative of the nonlinearity /. Thus the periodic waves 
depend on three parameters a, E, c together with a fourth constant of in- 
tegration Xo corresponding to the translation mode which can be modded 
out. Thus when we speak of a three parameter family of solutions we will be 
refering to a, E, c. On open sets in M 3 = (a, E, c) the solution to the above 
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is periodic. The boundary of this locus of points where the discriminant 
vanishes, which includes the constant solutions and the solitary waves. 

The constants a and c admit a variational interpretation: defining the 
Hamiltonian function 

H= / y - F{u)dx 
so that the Korteweg-DeVries equation can be written 

u t = ^5H(u) 

then the traveling waves are critical points of the augmented Hamiltonian 
functional 

5(H -aM - cP/2) = 

and thus represent critical points of the Hamiltonian under the constraint 
of fixed mass and momentum, with the quantities a, c representing Lagrange 
multipliers enforcing the constraints of fixed mass and momentum respec- 
tively 

We also note the connection with the underlying classical mechanics. 
The traveling wave ordinary differential equation is (after one integration) 
Hamiltonian and integrable. The classical action for this ordinary differential 
equation is 



K = j) pdq = y £ + aw + y F(u)du. 

The classical action provides a generating function for the conserved quan- 
tities of the traveling waves: specifically the classical action satisfies the 
following relationships 

-I 

In this paper, we are interested in both the spectral and orbital (non- 
linear) stability of spatially periodic traveling wave solutions of ([!]). The 
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spectral stability problems has been recently considered [HI EEU El] in which 
the authors considered stability to localized perturbations. In this case, the 
linearized stability takes the form 

JCfiv = jiv, (6) 

where C = —d x — f'(u) + c is a differential operator with periodic coefficients 
considered on the real Hilbert space L 2 (M) and J = d x gives the Hamiltonian 
structure. This is the standard form for the stability problem for solutions 
to equations with a Hamiltonian structure, although it must be emphasized 
that in the KdV case 3 = d x has a non-trivial kernel (spanned by 1) which 
complicates matters somewhat. 

Notation 1. Throughout this paper T will denote the period of the underlying 
periodic traveling wave. We will let 

T k = R/(fcTZ) 

be the torus of length kT and L 2 (Tk), the Hilbert space of square integrable 
kT periodic functions. We will frequently restrict to the subspace of mean 
zero functions. Following the notation of Deconinck and Kapitula we will 
denote this subspace by Hi : 

Ht = {0 G L\T k ) | (1,0) = 0}. 

In this paper we give a geometric construction of the spectrum of the 
linearized operator about a periodic traveling wave solution. This construction 
depends on the Jacobian determinants of various maps. Throughout this 
paper we will use the following Poisson bracket style notation for Jacobian 
determinants 

{F,G} Xj y = 

with the analogous notation for larger Jacobian determinants: 

F F F 

± x ± y ± z 

G x G y G z . 
H x H y H z 

We will also define the following eigenvalue counts: 



F F 

G Tj Gq, 



{F, G, H} Xi y tZ — 



1 INTRODUCTION AND PRELIMINARIES 



5 



Definition 1. Given the linearized operator d x C acting on Lp er (Tfc) we define 
ks, to be the number of real strictly positive eigenvalues, kc to be the number 
of complex-valued eigenvalues with strictly positive real part, and fcf to be 
the number of purely imaginary (non-zero) eigenvalues with negative Krein 
signature - in other words eigenvalues such that the corresponding eigenfunc- 
tions w satisfy 

(w,Cw) < 0. 

It is worth making a few remarks on this definition. The Krein signature is 
an important geometric quantity associated with eigenvalue problems having 
a Hamiltonian structure [37], and is associated with the sense of transversal- 
ity of the root of the eigenvalue relation. It is a fundamental result that 
two eigenvalues of like Krein signature collide they will remain on the axis, 
while if two eigenvalues of opposite Krein signature collide they will (gener- 
ically) leave the imaginary axis. It can be shown that a band of spectrum 
on the imaginary axis with multiplicity one will have eigenvalues of positive 
Krein signature, while a band of spectrum of multiplicity three will have 
two eigenvalues of positive Krein signature and one eigenvalue of negative 
Krein signature. It is not hard to check that outside of a sufficiently large 
ball in the spectral plane the spectrum lies on the imaginary axis and has 
multiplicity one: it follows from this that the non-imaginary eigenvalues and 
imaginary eigenvalues of negative Krein signature (being the roots of an an- 
alytic function) must be finite in number. Also notice that by symmetry 
of the spectrum about the real and imaginary axes, the quantities kc and 
kj are necessarily even, while the quantity &k has no definite parity. Our 
immediate goal is to establish the following index theorem, which is the main 
result of the paper: 

Theorem 1. Suppose u is a periodic solution to (JSJ) and T,M,P be the 
period, mass, and momentum of this solution considered as functions of the 
parameters (a, E, c) : 
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Also let fen, fee, k{ be the number of real, complex and imaginary eigenvalues 
of negative Krein signature of d x C on Lp er (Tfc) as defined above. Finally 
suppose that none of T E , {T, M} atE , {T, M, P} a ,E,c are zero. Then we have 
the following equality 



h +/cr+/cc = 2k— 1 




{TM,P} a 



E,c 



{T,M} aE 

{T,M,P}g. E ,C 

{TM} aE 
(10) 



The count on the left hand side of (flOj) clearly gives information con- 
cerning the spectral stability with perturbations in L^ er (T k ). In particular, 
if the count is zero then one can conclude the underlying periodic wave is 
spectrally stable to perturbations in Lp er (Tfc). Moreover, if the count is odd 
one is guaranteed the existence of a real eigenvalue since the quantities kc 
and kj are necessarily even. What is less clear is the relation of this count 
to the orbital stability of such a solution in L^ er {Yk). We will discuss this 
relationship in detail in the next section. The right hand side of (fTOl) relates 
to geometric information concerning the underlying periodic traveling wave 
profile, and the equality provides a direct relationship between the conserved 
quantities of the governing PDE flow restricted to the manifold of periodic 
traveling wave solutions to the stability of the underlying wave itself. 

In the case F is polynomial the integrals ([7]), ([8]), and ([9]) are Abelian 
integrals on a Riemann surface and the above expressions can be greatly 
simplified. For instance for the case of the Korteweg-de Vries equation 
the quantities Te, {T, M} a< E, {T, M, P} a ,E,c are homogeneous polynomials 
of degrees one, two and three resepctively in T, M, while for the modifed 
Korteweg-de Vries equation they are homogeneous polynomials in T and P. 
In general for a polynomial nonlinearity they are homogeneous polynomials 
of degree one, two and three in some finite number of moments of the solution 
Jfc = <£ -. == Thus, in the polynomial nonlinearity case Theorem 

[1] yields sufficient information for the stability of a periodic traveling wave 
solution in terms of a finite number of moments of the solution itself. 



< 
> 



2 Background and Main Results 

The study of eigenvalues of operators of the form ([6]) has a long history. 
Eigenvalue problems of exactly this form arise in the study of the stability 
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of solitary wave solutions to equations of Korteweg-de Vries type. The basic 
observation is that, if £ were positive definite the spectrum of 3C would 
necessarily be purely imaginary, since this operator is skew-adjoint under 
the modified inner product ((u,v)) = (C l l 2 u, £}l 2 v ), with (-, ■) the standard 
inner product on L 2 (IR). While in the case of nonlinear dispersive waves L is 
never positive definite due to the presence of symmetries one can count the 
number of eigenvalues off of the imaginary axis in terms of the dimensions 
of the kernel and the negative definite subspace of C. 

In the case of periodic solutions to the Korteweg-DeVries equation the 
best results of this type that we are aware of are due to Haragus and 
Kapitula[2T] and Deconinck and Kapitula[12j. Kapitula and Deconinck give 
the following construction: consider the spectral problem (ED with the lin- 
earized operator 3C acting on the real Hilbert space L 2 (T fc ), and let k^, kc, 
and ki be defined as before. Then one has the count 

K + k R + k c = n(C\ Hl ) - n(D). 

where n(-) denotes the dimension of the negative definite subspace of the 
appropriate operator acting on L 2 (Tfc), and D is a symmetric matrix whose 
entries are given by 

A,j = (yi,£\ Hl yj) 

where {yi} is any basis for the generalized eigenspace of d x C such that 

JC\ Hl span{?/i} = ker(£\ Hl ). 

The importance of this formula is the following: By using the results of |20j , 
it is known that a sufficient condition for the orbital stability of a periodic 
traveling wave solution of ([]]) is given by n(C\ Hl ) = n(D). Thus, if one is 
able to prove that 

K + k R + k c = 0, 

one can immediately conclude orbital stability in L 2 (l\.). It follows that 
spectral stability can be upgraded to the orbital stability if there are no 
purely imaginary eigenvalues of negative Krein signature. Moreover, the 
count clearly gives information concerning the spectral stability of the un- 
derlying periodic wave. In particular, a necessary condition for the spectral 
stability of such a solution in L 2 (Y k ) is for the difference n(£|^ 1 ) — n(D) to 
be even: more will be said on this later. 
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In another paper Bronski and Johnson [TT] considered the analogous spec- 
tral stability problem to localized perturbations from the point of view of 
Whitham modulation theory. Bronski and Johnson gave a normal form cal- 
culation for the spectral problem in a neighborhood of the origin in the 
spectral plane, which amounts to studying the spectral stability of a periodic 
traveling wave solution of ([T]) to long-wavelength perturbations: so called 
modulational instability. It was found that the presence of such an insta- 
bility could be detected by computing various Jacobians of maps from the 
conserved quantities of the gKdV flow to the parameter space (a, E, c) used 
to parameterize the periodic traveling waves. By deriving an asymptotic 
expansion of the periodic Evans function 

D(/j,e iK ) =det(M(/i)-e iK l) 

in a neighborhood of (//, k) = (0, 0), where M(/x) is the monodromy matrix 
associated with third order ODE (jSJ) and I is the three-by-three identity 
matrix, it was found that the spectrum of the operator JC in a neighborhood 
of the origin is determined by the modulational instability index 

Ami = \ ({T, P} E , C + 2{M, P} £ , a ) 3 - 3 Q{T, M, P} a , E )j ■ (11) 

In particular, it was found that if Ami > then the spectrum locally consists 
of a symmetric interval on the imaginary axis with multiplicity three, while 
if Ami < the spectrum locally consists of a symmetric interval of the 
imaginary axis with multiplicity one, along with two branches which, to 
leading order, bifurcate from the origin along straight lines with non-zero 
slope. Such information is important if one wishes to consider spectral or 
nonlinear stability to perturbations whose fundamental period is an integer 
multiple of that of the underlying wave. 

A similar geometric construction was later found useful by Johnson [23] 
to prove that such a solution with fundamental period T is orbitally stable 
in L 2 (Tk) if the quantities T E , {T, M} a ^ E , and {T, M, P} a ,E,c are all positive. 
Notice that for this sign pattern Theorem [1] also implies orbital stability in 
L 2 (Tfc), and hence the index theorem can be regarded as an extension of 
the theorem of Johnson. Both the calculations of Bronski and Johnson and 
of Johnson required a detailed understanding of the structure of the kernel 
and generalized kernels of the linear operators d x C and Cd x with periodic 
boundary conditions. These kernels were constructed by taking infinitesimal 
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variations in each of the defining parameters a, E, and c as well as the 
translation mode: this construction will be reviewed in the Proposition [TJ 
This information on the structure of the null-spaces, together with some 
additional work, will allow us to prove our main result. 

The operators we consider in this paper are non-self-adjoint and the null- 
spaces typically have a non-trivial Jordan structure. We adopt the following 
notation: 

Notation 2. Given an operator A acting on L 2 (Tk) for some k G N, we 
define the k th generalized kernel as follows 

g-ker fe (^) = ker{A k+1 )/ker{A k ). 

Thus g— ker (^4) = kei(A) is the usual kernel, A(g— ker 1 (A)) = ker(A), and 
so on. 

We begin by stating a preliminary lemma regarding the Jordan structure 
of d x C 

Proposition 1. Given any k G N, one generically has dim(ker(£)) = 1, 
dim(ker(d x C)) = 2, dim(g— ker 1 (9 x £)) = 1, and g— kexJd x C) = for j > 2. 
In particular, we have the following result: 

• IfT E ^ then ker£ = spanju^}. IfT E = then kerC = span-Ju^, Ue}- 

• If Te and T a do not simultaneously vanish then 



ker(d x £) = span < u x 



u E T e 

U n T n 



;i2) 



kei(Cd x ) = span{l,w} (13) 

// Te and T a simultaneously vanish then 

ker(d x £) = sp&n{u x , u a , u E } (14) 
kev(Cd x ) — span ^1, u, J ue dx^ . (15) 

Since the defining ordinary differential equation is third order the kernel 
cannot be more than three dimensional. 
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If{T,M} a , E ^0 then 



g— ker 1 ((9 x £) = span 
g— ker^Cdx) = span 



u E T e M e 
u a T a M a 
u c T c M c 

u E T e M e 



T a M a 

T r M r 



(16) 
(17) 



The generalized kernels g— ker 1 are one dimensional unless {T, M} atE 
and {T, P} a ,E vanish simultaneously. 

• Assuming {T, M} a E ^ the subsequent generalized kernels g— ker fc 
(k > 2) are empty as long as {T, M, P}e,o,,c 7^ 

Proof. This follows from the observation that the derivatives of the wave 
profile u with respect to the parameters a, E, c satisfy the following equations 



Cu x = 
Cue = 

Cu a = -1 (= 



Cu c = 



-u 



(= 



5M 
5u 

•-) 

5u ' 



) 



(18) 
(19) 

(20) 
(21) 



reflecting the fact that the constants (a, c) arise as Lagrange multipliers to 
enforce the mass and momentum constraints. In the above equality C denotes 
the formal operator without consideration for boundary conditions. In order 
to find elements of the kernel one must impose periodic boundary conditions. 
It is not hard to see that u x is periodic while derivatives with respect to the 
quantities are not periodic - since the period T depends on (E, a, c) "secular" 
terms (in the sense of multiple scale perturbation theory) arise: in particular 
one sees that the change across a period is proportional to derivatives of the 
period: 

/ u E {T) \ ( u E (0) \ ( u xE (0) \ 

u xE (0) 
u xxE (0) 



u xE (T) 
u xx e(T) 



\ 



J \ 



Tf 



J 



u xxE (0) 
u xxxE (0) 



\ 



/ 
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with similar expressions for the change in the u a ,u c across a period. Thus 
the quantity 

, / x ue Te 

<pi[x; a, c, e) = „ 

U>a J- a 

is periodic and satisfies C(f>i = Te- Similarly the quantity 

u E T e M e 



^ 2 (x;E,a,c) 



u a T a M a 
u c T c M c 



is by construction periodic and satisfies £0 2 = {T, M} E , C ~ {T, M} E , a u, and 
thus d x C(f)2 = {T, M}e,o,u x £ ker(d x C). Note that while 0i is essentially 
uniquely determined 02 is only determined up to an element of the kernel. 
Here we have chosen to make 02 have mean zero since this is the convention 
required in the work of Deconinck and Kapitula. More will be said on this 
choice later. 

The rest of the calculation follows in straightforward way from calcu- 
lations of this sort. For instance the existence of a second element of the 
generalized kernel is equivalent to the solvability of 



d x C = 



u E Ti 
T 

-L a. 



E 



U n 



By the Fredholm alternative ran((9 x £) = ker(£(9 x )- L and thus the above is 
solvable if only if 



(1, 

(u, 



u E 


T E 






u E 


T E 


U a 


T a 



= {T, M} a ,E = 

= {T,P} a , E = 



(22) 
(23) 



The rest of the claims follow similarly. In the case that the genericity con- 
ditions do not hold we do not attempt to compute the Jordan form, but we 
do remark that the algebraic multiplicity of the zero eigenvalue must jump 
from three to at least five, and is necessarily odd. □ 

In essence the above proposition shows that the elements of the kernel 
of d x C are given by elements of the tangent space to the (two-dimensional) 
manifold of solutions of fixed period at fixed wavespeed, while the element 
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of the first generalized kernel is given by a vector in the tangent space to 
the (three-dimensional) manifold of solutions of fixed period with no restric- 
tions on wavespeed. As one might expect all of the geometric information 
on independence in the above proposition can be expressed in terms of var- 
ious Jacobians. The next fact we note is that the signs of certain of these 
quantities conveys geometric information about the various operators. 

Lemma 1. Let n{C) be the dimension of the negative definite subspace of £ 
as an operator on L 2 (Yk) with periodic boundary conditions. Then 



Proof. As noted above u x G ker(£) and thus zero is a band-edge of C. From 
the Sturm oscillation theorem it is clear that either n(C) = n Ux — 1 (if zero 
is an upper band-edge) or n(C) = n Ux if zero is a lower band-edge. The 
Floquet discriminant k(/i) has positive slope at an upper band-edge and 
negative slope at a lower band-edge (and vanishes at a double point), and 
thus serves to distinguish the two cases. It can be shown (see appendix) that 
the sign of k'(fi) is equal to the sign of Tg and thus the result follows. Note 
that the sign of the slope of the Floquet discriminant has an interpretation 
as the Krein signature of the eigenvalues in the band. □ 

This lemma implies that the vanishing of Te signals an eigenvalue of C 
passing through the origin and a change in the dimension of n(C). The next 
proposition gives a similar interpretation for n(C\Hi)- 

Proposition 1. Let H\{= raxi{d x C) denote the space of mean-zero, T pe- 
riodic functions, and n(£|#J denote the dimension of the negative- definite 
subspace of the restriction of the operator C to this subspace. Assume that 
{T, M} a E 7^ and T E and{T, M} a ^ E never vanish simultaneously. Then we 
have the equality 




T E >0 
T e <0 




Proof. Since we are restricting to a codimension one subspace the Courant 
minimax principle immediately implies that we have either n(£|#j) = n(C) 
or n(C\ Hl ) = n(C) - 1. 
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Next note that when {T, M} a) E = then the function {u, T} aj E has mean 
zero and lies in H\. Further one has that 

£{u,T} atE = -T E 

or equivalently 

^\H 1 { U ^ T }a,E = 0. 

Therefore the vanishing of {T, M} Qj £ signals a change in the dimension of 
ti^CIh^- A local perturbation analysis shows that near a zero of {T, M} a< E 
we have that there is an eigenvalue of C\h x which is given by 

A = -^^k£ + 0({T ,M} a ,,) 

Thus as long as {T, M} a ^ vanishes for some parameter value the above 
count is correct. The case that {T, M} a} E is non-vanishing will be handled 
later. □ 

Remark 1. The last two results show that geometric quantities associated 
to the classical mechanics of the traveling waves contain information about 
changes in the nature of the spectrum of the linearized problem. Specifically: 

• Vanishing ofTs = Kee signals a change in the dimension ofn(C), the 
negative definite subspace of C. 



Vanishing of {T, M} a E 



K E e K Ea 
K aE K aa 

n{C\ Ran(9 x ))? the negative definite subspace of C-\~R&n(d x ) 



signals a change in the 



• Vanishing of {T, M, P} a ,E,c signals a change in the length of the Jordan 
chain of d x C. 

Finally, to conclude the proof of Theorem [H we must calculate n(D). 
This is the content of the following lemma. 

Lemma 1. Under the assumptions of Theorem^ one has that Del with 

D = -{T,M} a , E {T,M,P} a , E , c . 

Thus, n(D) is either or 1 depending if {T, M} ajE {T, M, P} a ,E,c is negative 
or positive, respectively. 



2 BACKGROUND AND MAIN RESULTS 



14 



Proof. Under the assumptions of Theorem[T], we know that ker(£) = span{ii x } 
and 

u E T e M e 



dM 



Hi 



U a 



T a 
T. 



M a 

M r 



{T, M} a}E u a 



from Proposition [IJ It follows that the matrix D is a real number in this 
case, with value equal to 



D 



u E 


T E 


M E 




u E 


T E 


M E 




Ta 


M a 


,£ 




T a 


M a 


U c 


T c 


M c 




u c 


T: 


M c 



■{T,M} a>E {T,M,P} a ,E,, 



as claimed. 

To complete the proof of proposition [I] note that the results of Bronski 
and Johnson imply the following identity: 



Omod 2 
lmod 2 



{T,M,P} a , E „ 
{T, M, P} a , E „ 



> 
< 



Since fcj" and kc are even they do not change the count modulo two. In the 
case {T, M} a E is non-vanishing the count is determined to within one, and 
is thus the count is exact if one knows the parity. Applying the result of 
Bronski and Johnson thus determines the count. □ 

Remark 2. It should be noted that the quantity D computed in Lemma 
also arose naturally in [23] when considering orbital stability of periodic trav- 
eling wave solutions of to perturbations with the same periodic structure. 
There, the negativity of D was necessary in order to prove the quadratic 
form induced by £ acting on L 2 {Tk) was positive definite on an appropriate 
subspace. As the methods therein are based on classical energy functional cal- 
culations, such a requirement was necessary to classify the periodic traveling 
wave as a local minimizer of the Hamiltonian subject to the momentum and 
mass constraints. 

The proof of Theorem [1] is now complete. Notice that Theorem [1] gives 
a sufficient requirement for a spatially periodic traveling wave of (OQ) to be 
orbitally stable in L 2 {T^) for any fceN and any sufficiently smooth nonlin- 
earity /. In the next two sections, we analyze Theorem [1] in the case of a 
power-nonlinearity by using complex analytic methods to reduce the expres- 
sion for the Jacobians involved in fflOl) in terms of moments of the underly- 
ing wave itself. This has the obvious advantage of being more susceptible 
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to numerical experiments as one no longer has to numerically differentiate 
approximate solutions with respect to the parameters a, E, and c. We will 
also discuss the computation of Ami for power-law nonlinearities. In partic- 
ular, we will prove a new theorem in the case of the focusing and defocusing 
MKdV which relates the modulational stability of a spatially periodic travel- 
ing wave to the number of distinct families of periodic solutions existing for 
the given parameter values. 

3 Polynomial Nonlinearities and the Picard- 
Fuchs System 

One major simplification of this theory occurs when the nonlinearity f(u) 
is polynomial. In this case the fundamental quantities (T, M, P) are given 
by Abelian integrals of the first, second or third kind on a Riemann surface. 
While we cannot give a detailed exposition of this theory here the basics are 
very straightforward. Suppose that P(u) = a +a 1 u+. . . a n u n is a polynomial 
of degree n. If the polynomial is of degree 2g + 2 or 2g + 1 then the quantity 
u J0L- is an Abelian differential on a Riemann surface of genus q. If we define 

the k th moment of the solution u(x) as follows: 



for any loop 7 in the correct homotopy class. For our purposes we are in- 
terested in branch cuts on the real axis though none of what will be said in 
this section assumes this. In the context of the stability problem one only 
needs I . . . J 4 , since T E , T a , . . . P c can all be expressed in terms of these five 
quantities, but the theory requires that one consider all such moments. The 
main observation is that the above integrals {Ik} are again Abelian integrals 
and thus can be expressed in terms of 

In practice the simplest way to do this is to use the identity 




then one obviously has 
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for m 6 (0..n — 1) and 



l P'{u)du 



2m 



m— 1 



j-tj+m-l 



2n/x m _i 



m— 1 



(24) 



(25) 



j=0 



for m G {0, 1, n}. This gives a linear system of 2n — 1 equations in 2n — 1 
unknowns 



2n-l. 
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The matrix which arises in the above linear systems is the Sylvester matrix 
of P{u) and P'{u). It is a standard result of commutative algebra that the 
Sylvester matrix of P(u) and Q(u) is singular if and only if the polynomials 
P and Q have a common root. In our case P{u) and P'(u) having a common 
root is equivalent to P(u) having a root of higher multiplicity In the case 
where P has a multiple root the a pair of branch points degenerate to a pole 
and the genus of the surface decreases by one. We will later work an example 
where this occurs. 

For a given polynomial it is rather straightforward to work these out, 
particularly with the aid of computer algebra systems. In this paper we 
did many of the more laborious calculations with Mathematica. [22] Some 
examples are presented in the next section. 

4 Examples 

4.1 The Korteweg -de Vries Equation (KdV-1) 

The Korteweg-Devries (KdV) equation 



u t + u xxx + [u 
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is, of course, completely integrable. The spectrum of the linearized flow 
can in principle be understood by the machinery of the inverse scattering 
transform, in particular by the construction via Baker-Akheizer functions 
detailed in the text of Belokolos, Bobenko, Enolskii, Its and Matveev[2|). 
Nevertheless this problem provides a good test for our methods, which we 
believe to be considerably simpler and easier to calculate than the algebro- 
geometric approach. 

Assuming a traveling wave u(x — ct) the ordinary differential equation 
integrates up to 



u: 



CU X ~\~ U xxx ~^ 

u xx = a + cu — u 
u 2 

- = E + au + c— 



u 2 ) x 



u 
~3 





-V'[ 

3 



U 



--: E- V(u) 



(26) 
(27) 

(28) 



where V(u) = V(u; a, c) denotes the effective potential for the Hamilto- 
nian system. A fundamental quantity is the discriminant of the polynomial 
R(u] a, E, c) := E — V(u; a, c), which is given by 

disc(i?(u; a, E, c)) = ^ (16a 3 + 3aV - 3QEac - 6Ec 3 - 36E 2 ) . 

Notice the KdV equation has periodic solutions if and only if disc(-R(w; a, E, c)) 
is positive. Moreover, by scaling (and possibly a map u \— > —u) the wave 
speed c can be assumed to be c = +1. 

In this case, the Picard-Fuchs system is the following set of five linear 
equations: 

1/3 \ 
c/2 1/3 

-1 
c -1 / 



( E 




a 




a 
E 
c 
a 




c/2 
a 

-1 

c 
a 



( h\ 
h 
h 
h 
h 



( T \ 

M 


2T 
\AM J 



where 



u 



du 



y/2R(u;a,E,c) 



and Ik 



u 



du 



(2R(u;a, E,c)) 



3/2 ' 



Solving this system implies the various Jacobians arising in Theorem [T] and 
the modulational stability index (fTTjl can be expressed in terms of the period 
T and the mass M as follows: 
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(4a + c 2 ) M + (6£ + ac) T 



12disc(i?(w, a, E, c)) 



{r,M} a , E 




{T, M, P} a ^ c 



12 disc(i?(w, a, E 1 , c)) 

r 3 (E-^(f)) 



2A M / 



2 disc(-R(w, a, £7, c)) 

(a 3 , T 3 + a 2 ,iT 2 M + a 1)2 TM 2 + a , 3 M 3 ) 2 



2 10 3 7 disc 3 (R(u, a, E,c)) 



where 



«0,3 



"3,0 



«2,1 



ai,2 



36£ + 18a£c - 8a 3 
18£c 2 - 6a 2 c + 36a£ 
(-18c£ + 24a 2 + 3ac 2 
c 3 + 6ac + 12£ 



These quantities are all positive. The positivity of Te follows from the result 
of Schaaf[3T] mentioned earlier. The non-negativity of Ami is clear: in prin- 
ciple the cubic polynomial in the numerator could vanish but numerics shows 
that it does not in the region where disc(_R(w, a, E, c)) > 0. The positivity of 
{T, M, P} a ,E,c is clear. Finally {T, M} a ^E is positive from Jensen's inequality 
since 



Remark 3. Recall from 01] / that the Jacobian {T, M, P} a ,E,c arises natu- 
rally as an orientation index for the gKdV linearized spectral problem for a 
sufficiently smooth nonlinearity. Indeed, one has that {T, M, P} a ,E,c < is 
sufficient to imply the existence of a non-zero real periodic eigenvalue of the 
linearized operator 3C, i.e. an unstable real eigenvalue in L 2 T [T/2, T/2]. 
Moreover, from fM^j it follows that if Te > 0, then such an eigenvalue can 
not exist if {T, M, P} a ,E,c is positive: however, no such claim can be made 
in the case where Te < 0. 

Theorem [1] now implies the following index result: if one considers the 
linearized operator acting on L 2 (Tfc) for k G N then u x has 2k roots in 
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Tfc and the number of real eigenvalues, complex eigenvalues, and imaginary 
eigenvalues of negative Krein signature satisfy 

fcj- + ku + k c = 2(k - 1) (29) 

In particular when k = 1, so one is considering stability to perturbations 
of the same period, the only eigenvalues lie on the imaginary axis and have 
positive Krein signature thus proving orbital stability of such solutions in 
L 2 (Tfc). Furthermore, considered as an operator on L 2 [— oo, oo] the spec- 
trum in a neighborhood of the origin in the spectral domain consists of the 
imaginary axis with multiplicity three, thus implying modulational stability 
of the periodic traveling wave solutions of the KdV equation. 



4.2 Example: Modified Korteweg- de Vries (KdV-2) 

The MKdV equation 

u t + u xxx ± {u 3 ) x = 

arises as a model for wave propagation in plasmas and as a model for the 
propagation of interfacial waves in a stratified medium. It is also integrable 
and the same caveats apply as for the KdV regarding the algebro-geometric 
construction of the spectrum of the linearized operator. The MKdV is invari- 
ant under the scaling x t— > ax,t \— > a 3 t,u i— > a"3«, and thus the wavespeed 
c can be scaled to be c = 0, ±1. The most physically and mathematically in- 
teresting case is the focusing MKdV (the plus sign above) with rightmoving 
waves where c can be scaled to +1. If we scale the focusing MKdV equation 
such that c = 1 the a, E parameter space contains the familiar swallow-tail 
fold: see Figure [H 

For the focusing MKdV, the swallowtail curve is defined implicitly by the 
equation 

disc(£ + au + u 2 /2 - w 4 / 4 ) = 2 « 2 ~ 27a 4 - 4E + 72a 2 E - 32E 2 - 64E 3 = 
or by the polynomial parametric representation 

a = s - s 3 (30) 

„2 o 4 

*-T"T- (31) 

The soliton solution corresponds to the origin a — 0, E — 0. Along the upper 
(dashed) branch (s G [— 3 , 3 ]) there are two solutions: a constant solution 
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Figure 1: The configuration space for focusing MKdV with c = +1. The 
swallowtail figure divides the plane into regions containing 0, 1, and 2 periodic 
solutions. The domain is colored according to the sign of {T, M, P} a ,E,c- 



and a solitary waves homoclinic to some (non-zero) constant value. Along the 
lower (dotted) branch s G [—1, — 3 ] U 1] there are again two solutions: a 
constant solution and a non-constant solution. Along the remaining portions 
of the curve there is only a constant solution. The Riemann surface associated 
to the traveling wave solutions is of genus one (a torus) 

ul „ u 2 u A 

— = E + au-\ 

2 2 4 

except along the swallowtail curve where the discriminant vanishes. In the 
case of vanishing discriminant the torus "pinches off" and degenerates to a 
cylinder. In this case all of the elliptic integrals can be evaluated in terms of 
elementary functions. Unlike the KdV case, where the only periodic solutions 
on the curve of vanishing discriminant are constant, the MKdV admits non- 
trivial periodic solutions on the swallowtail curve. 

The modulational instability index turns out, in this case, to be particu- 



4 EXAMPLES 



21 



larly simple. After solving the Picard-Fuchs system 
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one finds the following expressions for the various Jacobians: 

(3a 2 - IQE 2 - AE)T + (9a 2 - AE - l)M 



{T, M} a E 

{T,M,P} a ^ c 



16 disc(R(u; a, E, 1)) 
(3a 2 - AE)T 2 + (AE - 1)PT + P 2 

16 disc(i?(-u; a, E, 1)) 
(2a 2 - AE)T 3 + AEPT 2 - TP 2 + P 3 
32 disc(i?(w; a, P, 1)) 
(a , 3 T 3 + «i, 2 T 2 P + a 2jl TP 2 + aa.oi 33 ) 2 



where 



4194304 disc(P(w;a,P, l)) 3 
R{u; a, E, 1) := E + au + w 2 /2 - w 4 / 4 



as before and 



a 3 = (1 + 36P - 27a 2 ) 

a 2 ,i = (27a 2 + 144P 2 - 60P) 

a 1>2 = (36P - 240P 2 - 18a 2 + 108P) 

a 0)3 = (54a 4 - 180a 2 P + 144P 2 + 64P). 

We note a few things. First notice that, while the Picard-Fuchs system 
involves T, M, and P the resulting Jacobians only involve T and P. While 
this is not obvious from the point of view of linear algebra there is a clear 
complex analytic reason why this must be so: the Abelian differentials defin- 
ing T and P have zero residue about the point at infinity, as do Te, T a . . . P c , 
while M has a non- vanishing residue at infinity. Thus Tg, T a . . . P c must be 
expressible in terms of only T and P. 
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Secondly we note that while there are two distinct families of solutions 
inside the swallowtail they have the same orientation index {T, M, P} a ,E,c and 
modulational instability index Ami- This is special to MKdV, and follows 
from the fact that the integrals over one real cycle can be simply related to 
the integrals over the other real cycle via 



du f du 



7left V 1 H "> ^ 1 W J Tright 



Tie 



y/R(u: 


a, E, c) 


u 


du 


y/R(u: 


a, E, c) 


u 2 


du 




a, E, c) 



Tleft V U ! ^1 W ■'Tright 



y/R(u; 


a, E, c) 


u 


du 


^R(u: 


a, E, c) 


u 2 


du 




a, E, c) 



(32) 

+ 2V2n (33) 
(34) 



by deforming the contour onto the other cycle and picking up the contribution 
from the residue at infinity (which vanishes for T and P and is one for 
M). Since the orientation and modulational instability indices are built of 
derivatives of the above quantities these indices must be the same for both 
families of solutions. This observation extends the calculation of Haragus 
and Kapitula [21] for the zero amplitude waves to the periodic waves on the 
swallowtail curve. 

As in the KdV case the modulational instability index, which is a homo- 
geneous polynomial of degree 6 in T and P, can be expressed as the square of 
a homogeneous polynomial of degree 3 over an odd power of the discriminant 
of the polynomial R(u; a, E,c). A similar expression holds in the defocusing 
case, as well as for general values of c. The sign of this quantity is obviously 
the same as of the sign of the discriminant of the quartic, which is in turn 
positive if the quartic has no real roots or 4 real roots, and negative if the 
quartic has two only real roots. Thus we establish the following surprising 
fact: 

Theorem 2. The traveling wave solutions to the MKdV equation 

U t + U xxx : 

are modulationally unstable for a given set of parameter values if the polyno- 
mial 

E + au + cu 2 /2 ± u 4 /4 
has two real roots, and is modulationally stable if it has four real roots. 
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Remark 4. In the case of focusing MKdV, Theorem^ implies that if the pa- 
rameter values give rise to one periodic solution then this solution is unstable 
to perturbations of sufficiently long wavelength. If there are two periodic so- 
lutions then the spectrum of the linearization about one of these solutions in 
the neighborhood of the origin consists of the imaginary axis with multiplicity 
three. For the case of defocusing MKdV the situation is reversed: for a given 
set of parameter values there can be at most one periodic solution, which has 
no spectrum off of the imaginary axis in the neighborhood of the origin. 

Note that while this problem is in principle completely solvable using 
algebro-geometric techniques, Theorem [2] appears to be new. We suspect 
this is because the classical algebro-geometric calculations are sufficiently 
complicated that they are difficult to do in general. For examples of this sort 
of calculation see the original text of Belokolos et. al.jl] as well as the papers 
of Bottman and DeconinckjH] and Deconinck and Kapitula[12j. 

We now summarize the more interesting situation of the focusing MKdV 
in Figure [1] and below: 

(a) There are two families of solutions in this region. For both of these 
solutions the modulational instability index is positive and thus in a 
neighborhood of the origin the imaginary axis is in the spectrum with 
multiplicity three. Solutions in this region have Te > 0, {T, M} a> E > 0, 
and {T, M, P} a ,E,c > implying k^ + k R + k c = 2(k - 1). 

• The solutions in the remaining regions have a modulational instabil- 
ity index that is negative showing that they are always unstable to 
perturbations of sufficiently long wavelength. 

(b) In this region T E < 0, {T, M} a , E < 0, and {T, M, P} a>E , c < implying 
kj + ku + kc — 2k — 1. 

(c) In this region T E < 0, {T, M} a , E > 0, and {T, M, P} a ,E,c < implying 
kj +k^_+kc = 2k — 1. As one crosses between regions b and c the indices 
n(£)\ H i) and n(D) both increase (resp. decrease) by one, leaving the 
total count the same. 

(d) In this region T E < 0, {T, M} a , E > 0, and {T, M, P} a ,E,c > implying 
fcf + k u + k c = 2(k-l). 

(e) In this region T E > 0, {T, M} a , E > 0, and {T, M, P} a , E , c > implying 
kj+k R + k c = 2(k-l). 
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In regions (a), (d) and (e) when considering periodic perturbations (k — 1) 
one finds that fcf + /cr + kc = implying both spectral and orbital stability 
in L 2 (Ti). However, as mentioned above, in regions (d) and (e) the solution 
is spectrally unstable in L 2 (T^) for k e N sufficiently large. Moreover, in 
regions b and c there always exists a non-zero real periodic eigenvalue, i.e. the 
linearized operator 3C acting on L 2 (Tfc) always has a non-zero real eigenvalue 
and hence such solutions are always spectrally unstable. 

Remark 5. It should be noted that the above counts are consistent with the 
calculations of Deconinck and Kapitula IW{ in which they consider stability 
of the cnoidal wave solutions 

U(x, t; k) = \pl{i cn (yux — fi 2 (2 — K 2 )t; k) 

of the focusing MKdV equation, where /i > and k e [0, 1). Such solutions 
always correspond to regions (b) and (d) along with the constraint a = 0. 
There, the authors find numerically that there is a critical elliptic-modulus 
k* rs 0.909 such that solutions with < k < k* , corresponding to region (b) 
are orbitally stable in L 2 (Tf,) while solutions with k* < k < 1, corresponding 
to region (d) are spectrally unstable in L 2 (T k ) for allk e N due to the presence 
of a non-zero real eigenvalue of the linearized operator. 

4.3 Example: L 2 critical Korteweg-de Vries (KdV -4) 

Finally we consider the equation 

u t + u xxx + ty% = o. 

This equation is not a physical model for any system that we are aware of 
but is mathematically interesting for a number of reasons. This is the power 
where the solitary waves first go unstable. Equivalently this is the L 2 critical 
case, where one has the scaling u i— > ^u^x) preserving the L 2 norm and the 
relative contributions of the kinetic and potential energy to the Hamiltonian. 
Again we focus on the focusing case, which is the more interesting, and we 
scale everything so that c = +1. 

Again the parameter space is divided by a swallowtail curve into regions 
containing no periodic solution, one periodic solution, and two periodic so- 
lutions. The implicit representation is given by 

-48a 2 + 3125a 6 + 96£ - 11250a 4 £ + 10800a 2 £ 2 - 1728E 3 + 7776E 5 = 
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or parametrically by 



a = s — s 
_ s 2 _ 5s 6 
~ ~2 ~ IT' 



(35) 
(36) 



Again the picture is qualitatively similar to the MKdV case: the portion 
of the swallowtail parameterized by s 6 (— 5~*,5~4) represents parameter 
values for which there are two solutions: one constant and one homoclinic 
to a constant, with the origin representing the soliton solution (the solution 
homoclinic to zero) and the zero solution. The portions of the curves pa- 
rameterized by s G (— -^=,— l)U(i,l) represent parameter values for which 
there are two solutions: a constant and a periodic solution. The remainder 
of the curve represents parameter values for which there is only the constant 
solution. 

The Picard-Fuchs system is following set of eleven equations: 
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We have explicit expressions for the various Jacobians arising in the theory, 
but they are cumbersome and we will no reproduce them here: it should 
be noted they are homogeneous polynomials in /io(= T), M), /Z3, and 

/X4. Moreover, we note that these quantities turn out to be independent of 
P = /j,2 since the differential corresponding to momentum has a non-trivial 
residue at infinity, similar to the case of the MKdV. 

The stability diagram of these solutions is depicted in figure |2j Numerics 
indicate that the modulational instability index is always negative, indicating 
that solutions are always modulationally unstable. There are three curves 
emerging from the cusps of the swallowtail. The lowest of these is the curve 
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Figure 2: The configuration space for focusing KdV-4 with c = +1. The 
swallowtail figure divides the plane into regions containing 0,1, and 2 periodic 
solutions. 



on which the orientation index {T, M, P} a ,E,c vanishes, the middle (dashed) 
where T E vanishes and the upper (dotted) where {T, M} a ^ E vanishes. 
The behavior in the various regions is summarized as follows: 

(a) There are two solution families in this region, both of which satisfy 
T E > 0, {T, M} a , E > 0, and {T, M, P} a , E , c > 0. This implies that 
there are no real periodic eigenvalues and that kj + k^ + k c = 2{k — 1) 

(a') There is only one solution family in this region, otherwise the behavior 
is the same as in region (a') 

(b) The family of solutions in this region has T E > 0, {T, M} a ^ E > 0, and 
{T, M, P} a ,E,c < 0. This implies that /cf + k R + k c = 2k — 1 

(c) The family of solutions in this region has T E < 0, {T, M} a ^ E > 0, and 
{T, M, P} a ,E,c < 0. This implies that k^ + k R + k c = 2k — 1 . 

(d) The family of solutions in this region has T E < 0, {T, M} a>E < 0, and 
{T, M, P} a ,E,c < 0. This implies that fcf + k®. + k c = 2k — 1. 

It follows that solutions in region (a) are orbitally stable in L 2 (T^.) and 
spectrally unstable in L 2 (Tfc) for fceN sufficiently large. Moreover, solutions 
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in the remaining regions are spectrally unstable in L 2 (Tk) for any k G N due 
to the presence of a non-zero real periodic eigenvalue. 

5 Conclusions 

We have proven an index theorem for the linearization of Korteweg-de Vries 
type flows around a traveling wave solution and shown that the number of 
eigenvalues in the right half-plane plus the number of purely imaginary eigen- 
values of negative Krein signature given be expressed in terms of the Hessian 
of the classical action of the traveling wave ordinary differential equation or 
(equivalently) in terms of the Jacobian of the map from the Lagrange mul- 
tipliers to the conserved quantities. In the case of polynomial nonlinearity 
these quantities can be expressed in terms of homogeneous polynomials in 
Abelian integrals on a finite genus Riemann surface. 

The main drawback of the result is that it does not really distinguish 
between the eigenvalues in the right half-plane, which lead to an instabil- 
ity, and the imaginary eiegnvalues of negative Krein signature, which are 
generally not expected to lead to an instability. The index does distinguish 
between real eigenvalues and imaginary eigenvalues of negative Krein signa- 
ture, but only modulo two. This is sufficient to deal with the solitary wave 
case, where the only possible instability mechanism is the emergence of a sin- 
gle real eigenvalue from the origin. However in the periodic problem, where 
the behavior of the spectral problem is much richer, it would be preferable 
to have more information. 

We believe that a stronger result is possible: namely that there is spec- 
trum off of the imaginary axis if and only if the modulational instability index 
is negative. It follows from the results of this paper and the previous work 
of Bronski and Johnson that negativity of the modulational instability index 
is a sufficient condition for instability. Showing necessity would amount to 
showing that existence of any spectrum off of the imaginary axis would imply 
the existence of spectrum off of the imaginary axis in a neighborhood of the 
origin. In numerical experiments that we have conducted this has always 
been true. 
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